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ABSTRACT 

We have exploited ALMA calibration observations to carry out a novel, wide and deep submm 
survey, ALMACAL. These calibration data comprise a large number of observations of calibrator 
fields in a variety of frequency bands and array configurations. Gathering together data acquired 
during multiple visits to many ALMA calibrators, it is possible to reach noise levels which allow the 
detection of faint dusty, star-forming galaxies (DSFGs) over a significant area. In this paper we outline 
our survey strategy and report the first results. We have analysed data for 69 calibrators, reaching 
depths of ^ 25^Jybeam“^ at sub-arcsec resolution. Adopting a conservative approach based on 
> 5tT detections, we have found eight and 11 DSFGs in ALMA bands 6 and 7, respectively, with 
flux densities S'!. 2 mm > 0.2 mJy. The faintest galaxies would have been missed by even the deepest 
Herschel surveys. Our cumulative number counts have been determined independently at 870 fim 
and 1.2 mm, from a sparse sampling of the astronomical sky, and are thus relatively free of cosmic 
variance. The counts are lower than reported previously by a factor of at least 2 x. Future analyses 
will yield large, secure samples of DSFGs, with redshifts determined via detection of submm spectral 
lines. Uniquely, our strategy then allows morphological studies of very faint DSFGs - representative 
of more normal star-forming galaxies than conventional submm galaxies (SMGs) - in fields where 
self-calibration is feasible, yielding milliarcsecond spatial resolution. 

Subject headings: galaxy evolution; submm galaxies; dust emission; number counts 


1. INTRODUCTION 

Submm surveys revolutionised the study of galaxy 
formation and evolution by uncovering a population 
of dusty starbursts with submm flux densities of 
a few mJy, the so-called submm galaxies (SMGs , 
Smail et al.lll997t llvison et al.lll99^ iBarger et al.lll99§ : 
Hughes et al.l 1199'^ These relatively bright SMGs, 
which have a median redshift of z ^ 2.3 and SFR > 
lOOMQyr”^ (iChaoman et al.l 1200511 . have traditionally 
been found using bolometer cameras such as the Subm- 
millimetre Gqmmon -User Bolometer Array (SGUBA — 
iHolland et al.l 119^ mounted on single-dish telescopes 
such as the 15-m J a mes Clerk Maxwel l Telescope (e.g. 
Coppin et al.l l200^ iWeifi et al.l l2009bt IWardlow et al.l 

20111 iGeach et al.ll20I3l : iCasev et al.ll20I3[l . One of the 
main advantages of single-dish observations is that they 
can survey sufficiently wide areas of the sky to enable the 
detection of relatively large numbers of galaxies. How¬ 
ever, the large beam of single-dish telescopes makes it 
difficult to identify and explore the nature of even the 
brightest SMGs, unless dedicated interferometric follow¬ 
up observations are carried out, and any fainter popula¬ 
tion is buried in the confusion noise. 

Given that significant numbers of submm surveys have 
led to the selection of large samples of bright SMGs, 
the next obvious step is to carry out blind searches us¬ 
ing submm interferometers, detecting dusty star-forming 
galaxies (DSFGs) at sub-arcsec spatial resolution, at 
flux density levels below those accessible to single-dish 
telescopes. The main caveat to this approach is the 
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small field-of-view (FoV) of interferometric observations 
at wavelengths probing the dust emission from star¬ 
forming galaxies, such that even covering areas as small 
as a few arcmin^ is highly time intensive. 

To overcome this impediment, we are taking advantage 
of ALMA calibration scans to carry out a submm survey. 
Using calibrators it is possible to carry out a deep and 
wide submm survey with the necessary data coming ‘for 
free’ from science projects dedicated to a wide variety of 
astrophysical topics. A classical ALMA scheduling block 
(SB) comprises several steps, some of which involve ob¬ 
servations of very bright, compact sources with submm 
flux densities of the order of a Jy to calibrate the ampli¬ 
tude and phase of the visibilities of the science targets, to 
set the flux density scale and/or measure the bandpass 
response. Observations of such calibrators are essential 
and represent a significant fraction of each SB. A long list 
of calibrators is used to calibrate ALMA science datcQ. 
Each calibrator will typically be observed several times, 
on different dates, in several different ALMA bands, as 
part of one or several SBs corresponding to one or sev¬ 
eral ALMA science projects. By combining compatible 
data for a given calibrator we can reach r.m.s. noise lev¬ 
els sufficiently low to detect DSFGs within the primary 
beam (PB) centered on the calibrator. As an example, a 
typical observation of a bandpass calibrator lasts about 5 
minutes. According to the ALMA sensitivity calculator, 
in that time it is possible to reach a continuum depth of 
about 60 /iJy beam”^ in ALMA band 6 with 36 antennas, 
sufficient to detect DSFGs in the vicinity of the bandpass 

^ https://almascience.eso.org/sc/search 
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calibrator with flux densities, S'i. 2 mm > 0.3mJy at 5cr, 
ignoring for the moment the possible effect of limited 
dynamic range caused by the presence of a very bright 
source in the middle of the map. DSFGs are rare galax¬ 
ies, even at sub-mJy flux densities, so despite reaching 
very low noise levels, data for many calibrators must be 
acquired to increase the area surveyed and allow the de¬ 
tection of a significant number. 

There are several, key advantages of using calibrators 
to look for and analyse high-redshift DSFGs. ALMA 
calibrators tend to be observed in different projects with 
different ALMA configurations, ensuring excellent cover¬ 
age of the uv plane. Perhaps more importantly, a number 
of calibrators will be observed at extremely high spatial 
resolution, if this is amongst the requirements of the sci¬ 
ence project within which they are observed. The simul¬ 
taneous presence in the PB of one or more DSFGs and 
a bright ALMA calibrator lends itself perf ectly to self¬ 
calibration le.g. iPearson fc Readheadlll984 see also iQ]), 
which permits accurate tracking of the complex gains and 
hence near-perfect imaging, even with the longest avail¬ 
able interferometric baselines. This enables us to analyse 
the morphological properties of any fortuitously located 
DSFGs in unprecedented detail. The fact that each cal¬ 
ibrator is often observed in several different ALMA fre¬ 
quency bands allows us to ensure that faint detections 
are genuine, with close to 100% confidence (see 113.11) . 
Furthermore, such multi-band data allow us to study 
the spectral indices of DSFGs at matched spatial res¬ 
olution, and sometimes to determine their redshifts via 
so-called ‘blind’ detections of spectral lines such as CO. 
Next, the exposure time necessary for follow-up observa¬ 
tions of any DSFGs found using our approach are dra¬ 
matically reduced, since the calibrator field is used for 
the science observations, meaning that only a brief scan 
of a flux-density standard is required by way of calibra¬ 
tion. Finally, the number counts we report, coming from 
a sparse sampling of the astronomical s ky, are relatively 
free o f cosmic variance (see, for example. iSibthorpe et al.l 

[2(Ti^ . 

As part of our ALMA submm survey we aim to: (1): 
Study the submm number counts in different ALMA 
bands. So far, the number counts reported in previous 
work have been derived in different bands and conversion 
between flux densities has been required to provide a suf¬ 
ficient sample of galaxies. These conversion factors usu¬ 
ally assume the classical FIR/submm SED of an SMG at 
0 ^ 2.3 but not all submm detections are due to classical, 
dusty SMGs, and the redshift is rarely known for a faint 
source, meaning that these conversion factors between 
bands are uncertain, and consequently that multi-band 
information is the most robust way to derive number 
counts and carry out reliable comparisons with models 
of galaxy formation and evolution. (2): Search for emis¬ 
sion lines in the data cubes to determine the redshifts of 
the DSFGs via CO and other bright FIR/su bmm emis¬ 
sion lin es in the multi-band ALMA data ('e.g. lWeifi et al.l 
l2013al) . and to constrain the CO luminosity function and 
the cosmic H 2 de nsity, similar to the studies carried out 
so far with PdBI (jWalter et al.ll20l3: iDecarli et al.ll2014D 
but over much larger areas and using deeper observa¬ 
tions. (3): Carry out a morphological analysis of the 
emission from DSFGs in multiple bands at matched spa¬ 
tial resolution. 


In this work we focus on the description of the survey 
and the first derived number counts. Future work will 
describe the other aspects of the survey. 

This paper is structured as follows: in ^ we explain 
how the ALMA data used in this paper were obtained, 
calibrated and combined, ^details how those data have 
been analysed to construct our sample of DSFGs, includ¬ 
ing analyses of the contamination from calibrator-related 
sources, spurious sources, survey completeness, effective 
area covered and the effect of flux boosting. presents 
the number counts as determined at the current stage of 
the survey. Our conclusions are summarised in SjS] 

2. EXTRACTING AND CALIBRATING ALMA 
CALIBRATOR DATA 

In this section we explain how we retrieved and cal¬ 
ibrated the public ALMA calibrator data and how we 
created the maps we then used to detect submm emit¬ 
ters. 

We first used the ALMA calibrator archive to find all 
calibrators observed prior to 2015 July as part of science 
projects that passed quality control. Since we are pri¬ 
marily interested in the selection of star-forming galaxies 
via the detection of their redshifted dust emission, we re¬ 
trieved data in ALMA bands 6 (B6, around 250 GHz or 
1.2mm) and 7 (B7, around 345 GHz or 870/rm), where 
we then select sources (see 113.11) . The projects contribut¬ 
ing to our sample represent a random selection of all the 
projects undertaken at ALMA, so the biases in terms 
of sky coverage are generated only by the latitude of 
ALMA, the annual weather patterns in the Atacama 
desert, and the positions of objects of interest to the 
astronomical community, e.g. the Milky Way. 

For datasets outside the proprietary time period, we 
retrieved the full data deliveries from the ALMA archive. 
For datasets that remained within the proprietary time, 
a ticket to the ALMA helpdesk was submitted with a 
request to obtain the part of the data that includes the 
calibrator scans. 

During the execution of the ALMA calibration scripts, 
the bandpass calibrators are not always fully calibrated, 
in the sense that the calibration tables obtained from 
the so-called phase calibrators are always applied to the 
science targets but not always to the bandpass calibra¬ 
tors. For this reason, we needed to recover the correct 
flux density scale of the bandpass calibrators from the 
flux tables in the data delivery packages and re-apply 
the calibration tables to the calibrators. 

Next, we created so-called ‘pseudo-continuum’ mea¬ 
surement sets for which all channels in each spectral win¬ 
dow were averaged. These pseudo-continuum hies were 
used to self-calibrate the calibrator data. Two rounds of 
self-calibration were applied, at hrst only in phase, then 
in both amplitude and phase, both with a solution in¬ 
terval equal to the integration time. Instead of imaging 
the data inbetween the self-calibration steps, we used a 
point-source model that we htted to the uv data. The 
advantage of subtracting the point-source model for each 
observation separately is that any variability of the cali¬ 
bration source will not affect the calibration of the com¬ 
bined data. For the majority of the datasets we found 
that these two rounds of self-calibration and uu-model 
htting produced adequate results. Finally, the point- 
source model was subtracted from the visibility data. 















ALMACAL I: a deep and wide submm survey 


3 


This procedure produced calibrated visibilities for the 
background region of each calibrator scan. 

The calibrator-subtracted visibilities for every SB were 
then imaged individually, without combining data for a 
given calibrator, using the CLEAN task. To do this, we 
defined a cleaning window of radius 1.5 x the expected 
FWHM of the PB (23" and 17" in B6 and B7, respectively) 
in each band and cleaned down to the r.m.s. of the dirty 
image. We then inspected every map by eye, discarding 
all datasets which showed evidence of poor calibration. 
We also inspected the calibrated visibilities to discard 
any poorly-calibrated SBs. This led to the loss of a sig¬ 
nificant number of datasets (~ 20%), but their size and 
complexity made flagging and re-calibration impractica¬ 
ble. 

Having discarded all bad datasets, we re-calculated 
the weights of the visibilities of the remaining datasets 
using STATWT to measure the visibility scatter empiri¬ 
cally as a function of time, antenna, and/or baseline for 
each dataset, re-weighting the visibilities accordingly, be¬ 
fore combining all data for each calibrator in each band 
with CONCAT. We then created the deep maps follow¬ 
ing the same cleaning technique explained earlier. Ta¬ 
ble [T] provides a summary of the calibrators used in this 
survey, along with the depth and beam size of each in¬ 
dividual map. We have determined the FWHM of each 
calibrator map by fitting a Gaussian profile to the PB 
response obtained during the cleaning process in CASA. 
We have used 69 pointings, covering a total area close to 
19 arcmin^ for bright DSFGs, S'sTO/xm > 1 mJy, although 
the effective area for a given flux density depends on the 
depth of the maps - see ^13.71 It can be seen in Table [T] 
that most of our ALMA maps have sub-arcsec resolution 
and reach r.m.s. noise levels, a ~ 30/rJybeam“^, in B6 
and B7. 

We anticipated that one possible limitation of our ap¬ 
proach might be the presence of the bright calibrator in 
the middle of the map, which may have been expected to 
influence the dynamic range of the image, defined as the 
ratio between the flux density of the brightest source de¬ 
tected and the r.m.s. of the map. However, the r.m.s. of 
the clean maps listed in Table [T] indicates that we have 
reached dynamic ranges in excess of 18,000, detecting 
DSFGs as faint as 5'i.2 mm 0.2 mJy at the present stage 

of the survey (see 113.21) . We conclude, therefore, that the 
presence of the calibrator in the middle of the image has 
not proved to be a strong limitation. 

3. ANALYSIS 
3.1. Source detection 

As has been the case with most previous reports 
of faint ALMA counts le.g. I^tsulpde et al.l 120131 
Ono et ail 12014 iFuiimoto et al.l 120151 iCarniani et al l 
20151) . our source detection has be en performed using 
SExtractor (|Bertin &: Arnoutill996D on the clean maps, 
before correcting for the antenuation response of the 
PB. For those calibrators for which we have both B6 
and B7 data, we searched for sources in B6 because the 
maps have a larger FoV. Searching for sources in B7 pro¬ 
vided no additional sources. Because calibrators are such 
bright sources, we might expect some low-level residual 
emission in the maps, even after subtracting the calibra¬ 
tors from the data using a point-source model. Faced 


with this fear, we adopted a conservative source detec¬ 
tion technique to exclude the possibility of spurious de¬ 
tections: we selected sources with peak flux densities at 
least 5x the r.m.s. noise, rather higher than the SNR ~ 4 
limit used in previous work. Our higher SNR threshold 
means we do not detect very faint sources, but we reduce 
the number of false detections considerably. Despite our 
conservative approach, we were able to detect sources 
down to SWoum ^ 0.4 mJy at sub-arcsec spatial resolu¬ 
tion (see )13.2p . even at the present stage of the survey. 

After correcting for the PB response, the r.m.s. of the 
clean maps increases rapidly with the distance from the 
map center and hence our sensitivity to sources at large 
distances from the map center decreases. We searched 
for detections within a radius of 1.5 x the fwhm of the 
PB of each map, a region which is not severely affected 
by PB attenuation. A « 1-mJy galaxy was detected close 
to the edge of the detection area around the calibrator, 
J1744—3116, showing that this detection area does not 
include inaccessible parts of the PB-corrected maps. 

Whenever we detected an DSFG around a calibrator, 
the ALMA archive was queried again, this time with no 
limitation placed on the ALMA observing band. This 
additional search is done to: (1) increase the SNR of 
the detection; ( 2 ) confirm via multi-band observations 
that faint sources are real; (3) look for emission lines in 
the FIR/submm spectrum of the submm emitter to de¬ 
termine or constrain its redshift; (4) exploit multi-band 
information to distinguish between genuine DSFGs and 
jets emanating from the calibrator or other calibrator- 
related emission. Data obtained in ALMA bands 3 and 
4 ofte n allow the detec t ion of the brightest mid-J CO 
lines (I Weiss et al.l l2007t iWeifi et all l2009aL l2013bll and 
also represent the most efficient way to distinguish be¬ 
tween synchrotron-powered jets and thermal emission 
from DSFGs; data obtained in ALMA bands 8 and 9 can 
improve our sampling of the FIR/submm SED - at least 
for sources close to the calibrator, within the smaller, 
high-frequency EoV. 

3.2. Source catalogue 

In our 69 ALMA maps we have found eight and 11 
submm detections in B6 and B7, respectively. Account¬ 
ing for all the galaxies detected in B6 and B7, this rep¬ 
resents a total sample of 13 submm detections at > 5a. 
The ten calibrator fields in which our submm emitters 
were found are shown in Figure [TJ 

To ensure that the flux densities of the detected sources 
are well determined, we re-imaged the visibilities of the 
calibrators around which our DSFGs are located using 
a slightly different procedure. We defined around each 
submm detection a cleaning box with 1.5" on each side, 
a value greater than all our synthesised beams (see Ta¬ 
ble [T]) and cleaned down to the r.m.s. of each dirty map. 
As will be explained in 113.81 this non-interactive clean¬ 
ing method provides the most accurate determination of 
the flux density of the detected DSFGs. After cleaning, 
the maps were corrected for PB attenuation using imp- 
BCOR and the flux densities and uncertainties were then 
determined in the PB-corrected maps using imfit, with 
the same box used during cleaning. The coordinates and 
multi-band flux densitieso of the detected DSFGs are pre- 

^ Quoted uncertainties include the fitting errors; since we do 
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TABLE 1 

Summary of the observations at the present stage of the survey 


Calibrator 

R.A. 

[J2000] 

Dec. 

[J2000] 

z 

[^Jy beam“^] 

,-, 

'—' 

[/i,Jy beani”^] 

["] 

JOOll-2612 

00:11:01.2 

-26:12:33.1 

1.096 

47 

0.35 

X 

0.30 





J0121+1149 

01:21:41.6 

+11:49:50.4 

0.570 

85 

0.52 

X 

0.50 





J0132-1654 

01:32:43.5 

-16:54:48.5 

1.020 

84 

0.45 

X 

0.34 





J0141-0928 

01:41:25.8 

-09:28:43.7 

0.733 

98 

0.37 

X 

0.22 





J0145-2733 

01:45:03.4 

-27:33:34.3 

1.155 

91 

0.28 

X 

0.23 





J0238+1636 

02:38:38.9 

+16:36:60.1 

0.940 

47 

0.41 

X 

0.35 





J0239-0234 

02:39:45.5 

-02:34:41.0 

1.116 

75 

0.55 

X 

0.51 





J0329-2357 

03:29:54.1 

-23:57:08.8 

0.895 

73 

0.31 

X 

0.22 





J0519-4546 

05:19:49.7 

-45:46:43.9 

0.035 

39 

0.33 

X 

0.24 





J0550-5732 

05:50:09.6 

-57:32:24.4 

2.001 

93 

0.37 

X 

0.24 





J0854+2006 

08:54:48.9 

+20:06:30.6 

0.306 

53 

0.73 

X 

0.69 





J1215+1654 

12:15:04.0 

+16:54:38.0 

1.132 

92 

1.43 

X 

0.93 





J1308-6707 

13:08:17.2 

-67:07:05.0 


93 

0.33 

X 

0.27 





J1342-2900 

13:42:15.3 

-29:00:41.8 

1.442 

87 

0.68 

X 

0.51 





J1550+0527 

15:50:35.3 

+05:27:10.5 

1.422 

39 

0.63 

X 

0.45 





J1826-2924 

18:26:20.6 

-29:24:25.0 


67 

0.53 

X 

0.43 





J1832-1035 

18:32:20.8 

-10:35:11.2 


37 

0.50 

X 

0.37 





J1933-6942 

19:33:31.2 

-69:42:58.9 

1.481 

96 

0.85 

X 

0.50 





J1955+1358 

19:55:11.6 

+13:58:16.2 

0.743 

88 

0.61 

X 

0.51 





J2009-4849 

20:09:25.4 

-48:49:53.7 

0.071 

49 

0.52 

X 

0.43 





J2223-3137 

22:23:21.6 

-31:37:02.1 


37 

0.59 

X 

0.52 





J2306-0459 

23:06:15.3 

-04:59:48.3 

1.139 

67 

0.50 

X 

0.36 





J2331-1556 

23:31:38.6 

-15:56:57.2 

1.153 

32 

0.43 

X 

0.31 





J2357-5311 

23:57:53.2 

-53:11:14.0 

1.006 

78 

0.64 

X 

0.39 





J0217+0144 

02:17:49.0 

+01:44:49.7 

1.715 





44 

0.48 

X 

0.34 

J0339-0146 

03:39:30.9 

-01:46:35.8 

0.852 





92 

0.35 

X 

0.30 

J0607-0834 

06:07:59.7 

-08:34:50.0 

0.872 





41 

0.34 

X 

0.30 

J1048-1909 

10:48:06.6 

-19:09:35.7 

0.595 





65 

1.27 

X 

0.58 

J1303-5540 

13:03:49.2 

-55:40:31.6 






57 

0.37 

X 

0.29 

J1347+1217 

13:47:33.4 

+ 12:17:24.2 

0.122 





40 

0.61 

X 

0.42 

J1505+0326 

15:05:06.5 

+03:26:30.8 

0.408 





93 

0.46 

X 

0.29 

J1625-2527 

16:25:46.9 

-25:27:38.3 

0.786 





74 

0.52 

X 

0.29 

J1700-2610 

17:00:53.2 

-26:10:51.7 






97 

0.64 

X 

0.41 

J1744-3116 

17:44:23.6 

-31:16:36.3 






60 

0.33 

X 

0.29 

J1924-2914 

19:24:51.1 

-29:14:30.1 

0.353 





30 

0.33 

X 

0.29 

J2206-0031 

22:06:43.3 

-00:31:02.5 

0.335 





48 

0.36 

X 

0.31 

J2232+1143 

22:32:36.4 

+ 11:43:50.9 

1.037 





41 

0.41 

X 

0.37 

J0038-2459 

00:38:14.7 

-24:59:02.5 

0.498 

86 

0.49 

X 

0.39 

231 

0.43 

X 

0.29 

J0108+0135 

01:08:38.8 

+01:35:00.8 

2.099 

62 

0.69 

X 

0.49 

129 

1.02 

X 

0.56 

J0215-0222 

02:15:42.0 

-02:22:56.8 

1.178 

58 

0.65 

X 

0.52 

32 

0.39 

X 

0.34 

J0217+0144 

02:17:48.9 

+01:44:50.0 

1.715 

52 

0.31 

X 

0.25 

54 

0.33 

X 

0.27 

J0224+0659 

02:24:28.4 

+06:59:23.3 

0.511 

75 

0.72 

X 

0.59 

114 

0.36 

X 

0.31 

J0241-0815 

02:41:04.8 

-08:15:20.8 

0.005 

28 

0.54 

X 

0.37 

36 

0.41 

X 

0.36 

J0334-4008 

03:34:13.7 

-40:08:25.4 

1.445 

35 

0.34 

X 

0.25 

72 

0.29 

X 

0.21 

J0348-2749 

03:48:38.1 

-27:49:13.6 

0.991 

85 

0.37 

X 

0.23 

59 

0.28 

X 

0.22 

J0423-0120 

04:23:15.8 

-01:20:33.1 

0.916 

29 

0.46 

X 

0.44 

61 

0.33 

X 

0.23 

J0510+1800 

05:10:02.4 

+18:00:41.6 

0.416 

33 

0.54 

X 

0.42 

65 

0.32 

X 

0.23 

J0522-3627 

05:22:58.0 

-36:27:31.0 

0.057 

60 

0.74 

X 

0.55 

123 

0.44 

X 

0.28 

J0538-4405 

05:38:50.3 

-44:05:08.9 

0.890 

46 

0.42 

X 

0.32 

62 

0.61 

X 

0.46 

J0635-7516 

06:35:46.5 

-75:16:16.8 

0.653 

23 

0.37 

X 

0.28 

48 

0.32 

X 

0.21 

J0825+0309 

08:25:50.3 

+03:09:24.5 

0.506 

39 

0.50 

X 

0.42 

56 

1.05 

X 

0.58 

J0909+0121 

09:09:10.1 

+01:21:35.6 

1.025 

68 

1.28 

X 

0.89 

48 

1.14 

X 

0.55 

J1008+0621 

10:08:00.8 

+06:21:21.2 

1.720 

30 

0.50 

X 

0.48 

51 

0.28 

X 

0.23 

JlOlO-0200 

10:10:51.7 

-02:00:19.6 

0.890 

29 

1.25 

X 

0.73 

39 

0.91 

X 

0.51 

J1037-2934 

10:37:16.1 

-29:34:02.8 

0.312 

45 

1.30 

X 

0.70 

93 

0.80 

X 

0.50 

J1058+0133 

10:58:29.6 

+01:33:58.8 

0.890 

44 

0.50 

X 

0.47 

82 

0.31 

X 

0.23 

J1215-1731 

12:15:46.8 

-17:31:45.4 


83 

0.94 

X 

0.43 

54 

0.38 

X 

0.28 

J1229+0203 

12:29:06.7 

+02:03:08.6 

0.158 

171 

0.90 

X 

0.68 

80 

0.62 

X 

0.43 

J1337-1257 

13:37:39.8 

-12:57:24.7 

0.539 

83 

0.71 

X 

0.52 

50 

0.56 

X 

0.43 

J1427-4206 

14:27:56.3 

-42:06:19.4 

1.522 

35 

0.57 

X 

0.47 

43 

0.38 

X 

0.32 

J1517-2422 

15:17:41.8 

-24:22:19.5 

0.049 

93 

0.69 

X 

0.30 

42 

0.38 

X 

0.32 

J1534-3526 

15:34:54.7 

-35:26:23.0 

1.515 

50 

0.51 

X 

0.43 

116 

0.46 

X 

0.35 

J1617-5848 

16:17:17.9 

-58:48:07.9 


51 

0.38 

X 

0.24 

154 

0.46 

X 

0.37 

J1733-1304 

17:33:02.7 

-13:04:50.0 

0.902 

45 

1.06 

X 

0.55 

42 

0.39 

X 

0.36 

J1832-2039 

18:32:11.0 

-20:39:48.2 

0.103 

87 

0.43 

X 

0.38 

50 

0.40 

X 

0.33 

J2056-4714 

20:56:16.3 

-47:14:48.5 

1.489 

26 

0.47 

X 

0.44 

41 

0.56 

X 

0.47 

J2148+0657 

21:48:05.5 

+06:57:38.6 

0.791 

34 

0.51 

X 

0.46 

56 

0.31 

X 

0.23 

CTA102 

22:32:36.4 

+11:43:50.8 

1.037 

35 

0.55 

X 

0.42 

34 

0.55 

X 

0.42 

J2258-2758 

22:58:06.0 

-27:58:21.2 

0.926 

25 

0.36 

X 

0.30 

94 

0.38 

X 

0.31 
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Fig. 1.— ALMA images of the 13 DSFGs detected around calibrators (represented by the black squares). The calibrators lie at the 
center of each map, represented by the red cross; they have been subtracted in the uv plane, using point-source models, prior to imaging. 
These are 870-pm (ALMA band-7) images, except in the two cases where only band-6 data are available, shown prior to the correction for 
PB attenuation. For J0108+0135 we also use the ALMA band-6 image in order to show ALMACAL J010838.75-1-013455.9, part of the jet 
emanating from the calibrator as revealed by ALMA band-3 imaging. Each image is 25'' on each side (~ 1.5x the fwhm of the band-7 
PB). The jet emanating from the calibrator, J2223—3137, is clearly visible. N is up; E is to the left. 
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seated in Tabled] The S'a.omm column shows the ALMA 
B3 coverage of our sample, where available. None are 
detected (expressed by the symbol < cr), meaning their 
SEDs are all consistent with DSFGs. 

Most of our submm emitters are fainter in B7 (i.e. 
at around 870/rm) than typical SCUBA- or LABOCA- 
detected SMGs. Ten of our 13 DSFGs are fainter than 
the faint est deboosted SMG reported in the LESS survey, 
3.5 mJy (iWeifi et al.ll2009bL where the detection thresh¬ 
old was SNR > 3.75). Around half of our DSFGs are 
fainter than 1 mJy at 870 fim and they are all detected 
at relatively high SNR and with sub-arcsec spatial res¬ 
olution. The median 870-/im flux density of our DS¬ 
FGs is comparable with th e flux density of Ly man-break 
galaxies (LBGs) at z ^ 3 (|CoDPin et al.ll201^ . detected 
only via stacking. Since we aim to carry out a deep 
and wide survey, the faintest galaxies detected are of 
special interest. Figure [D shows the FIR SEDs of the 
faintest galaxies detected in B6 (left panel) and B7 (right 
panel), respectively. We do not yet know their redshifts, 
so we show SEDs for z = 0-6, adopti ng the shape of 
the median SED of the LESS SMGs (|Swinbank et al.l 
120141 . Figure [2| indicates, then, that our faintest DS¬ 
FGs are fainter than FIR-detected LBGs, F IR-detected 
sBzK galaxies or FIR-detected Ha emitters (|Oteo et al.l 
IMl EoH) at z < 2. Indeed, our faintest DSFGs 
would not have been detected by the deepest survey car¬ 
ried out by Herschel and they may thus represent the link 
between the extreme population of Herschel or SCUBA- 
/LABOCA-selected galaxies and the less extreme UV- 
selected galaxies. It is clear that our survey represents a 
significant step towards the discovery and characterisa¬ 
tion of faint, sub-mJy DSFGs. 

There are also a number of bright DSFGs in our sam¬ 
ple, with Ssrotim > 4mJy. The emission from these 
galaxies is confined to ~0.3", being only slightly re¬ 
solved, as suggested by their peak to integrated flux 
densities. They represent a population of extremely IR- 
bright galaxies whose prodigious star formation must 
be confined to a remarkably small vol ume (see also 
iSimpson et al.l[^15at llkarashi et al.l[2?)T^ . 

Among the 13 detected DSFGs, six have measure¬ 
ments in both B6 and B7. Figure [3| shows that their 
5'870/jm/<S'i.2mm ratios are generally in good agreement 
with those expected for DSFGs, as represented by the 
average S ED of SMGs in the LE SS survey at redshifts, 
z = 1-3 (|Swinbank et al.l I20l4l . One of the galaxies 
has an elevated 5'870/jm/<S'i.2mm which can be explained 
with an emissivity (3 > 2, already found for some DS¬ 
FGs (see for example iGasev et al.l 1201111 . The lowest 
'S' 870 /jm/<S'i. 2 mm ratio in our sample is commensurate with 
that galaxy lying at z > 3. This is the first time that 
FIR/submm spectral indices have been derived for bright 
or faint DSFGs at matched, high spatial resolution. The 
upper limits in 5 'b7,b6 /5'b 3 are also compatible with the 
SEDs of high-redshift DSFGs. 

Despite the smaller area covered by our B7 observa¬ 
tions (Figure |4|) we have detected more galaxies in B7 
than in B6. For the fields where data are available only 
in a single band, we still detect more galaxies in B7 than 
in B6. This apparent discrepancy is due to a combination 
of different factors, including the different noise levels in 
the maps, the low-number statistics at the present stage 
of the survey, the selection function of faint DSFGs as 


a function of wavelength, and the multiplicity of DSFGs 
( 113.31) . Among the five DSFGs detected only in B7, four 
are found in pairs, so the five galaxies are actually de¬ 
tected in only three maps. This is still larger than the 
number of maps with B6-only DSFGs (two), and merely 
reflects the need to survey large areas down to low noise 
levels to derive robust number counts. 

3.3. Multiplicity 

Previous studies exploring the environments of bright 
SMGs have reported that they tend to be strongly clus¬ 
tered (|Blain et al.l 12004 iScott et ahl l2006t iHickox et al.l 
l2012f) . Furthermore, it has been reported that 
bright SMGs found with single-dish telescopes are of¬ 
ten resolved into several different components when 

they ar e observed at _ hig h spatial resolution with 

ALM A (iKarim et al.l 120131 ISimpson et al.l I2015bll : in¬ 
deed, llvison et all (1200711 saw the same effect with the 
Very Large Array, arguing that many of the pairs must 
be physically associated rather than chance alignments 
of dusty galaxies in the same line of sight. 

Among our 13 DSFGs, six are found in pairs, with typ¬ 
ical separations of a few arcsecs, consistent with these 
earlier findings. Observed with single-dish telescopes, 
these pairs would be seen as single, unresolved emis¬ 
sion, similar to the classical population of SGUBA- or 
LABOGA-selected SMGs. At the present stage of our 
survey, we can confirm only that the two DSFGs seen 
around J1058-I-0133 are at the same redshift (Oteo et al. 
2015, in preparation). One of the advantages of using 
ALMA calibrators to study the DSFG population is that 
the calibrators will eventually be observed in a frequency 
range where there is one or more emission lines. A single 
line will allow us to confirm whet her the multiple compo - 
nents lie at the same redshift fe.g. lGhapman et al.ll201^ . 

3.4. Caveats 

As we have described, there are a number of advantages 
relating to the use of ALMA calibration data to study 
DSFGs. Nevertheless, we must note some caveats. 

The first relates to the possibility that unresolved, 
high-redshift DSFGs may be confused with jets emanat¬ 
ing from the calibrators. ALMA calibrators are typically 
blazars with strong jets that are currently oriented along 
our line of sight. Indeed, jet signatures are clearly visible 
in some of our ALMA maps. Jets tend to present an ex¬ 
tended, cometary shape and sometimes curved tails ema¬ 
nating from the calibrator (see, for example, J2223—3137 
in Figure [IJ. Due to their morphology, these jets are rel¬ 
atively easy to identify and thus distinguish from high- 
redshift DSFGs by visual inspection of the clean maps, 
since their extended emission points towards the calibra¬ 
tor. However, there may also be cases where the jet com¬ 
ing from a calibrator may appear unresolved. This un¬ 
resolved emission could be confused with a high-redshift 
DSFG. It is, however, possible to discriminate between 
jets and DSFGs using the spectral index of the emission 
in the ALMA data. The ratio between ALMA B6 or B7 
and B3 or B4 flux densities provide useful diagnostics to 
distinguish between thermal and non-thermal emission: 
'S'b 3,4 > *S'b 6,7 and S'b 7 < Sbg would be typical for a jet. 
When there are no ALMA B3 or B4 data, the flux den¬ 
sity ratio between B6 and B7 can also be used to identify 
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TABLE 2 

DSFGS DETECTED UP TO 2015 JULY IN OUR ALMA SUBMM SURVEY (SEE FIGURE [!}. 


Source 

Ssroiim [mJy] 

Si.2mm [mJy] 

Sa.Omm [mJy] 

ALMACAL J130349.05- 554034.2 

0.40 ±0.09 



ALMACAL J130349.43-554028.5 

0.66 ±0.08 



ALMACAL J174422.69-311639.4 

1.12 ±0.27 


< (T 

ALMACAL J220642.87-003108.1 

6.89 ±0.36 



ALMACAL J220642.98-003110.8 
ALMACAL J 193329.46-694258.4 

0.71 ±0.15 

1.42 ±0.20 


ALMACAL J222321.73-313707.7 


0.23 ±0.04 


ALMACAL J101051.03-020018.8 

0.52 ±0.15 

0.20 ±0.04 


ALMACAL J214806.00-1-065736.2 

2.09 ±0.32 

0.55 ±0.07 

< O 

ALMACAL J024104.82-081515.0 

0.89 ±0.10 

0.49 ± 0.06 

< O 

ALMACAL J010838.56-1-013504.2 

2.20 ±0.15 

0.80 ±0.17 

< (7 

ALMACAL J105829.54-I-013359.7 

6.48 ± 0.30 

2.16 ± 0.17 

< (7 

ALMACAL J105829.73-I-013357.2 

4.35 ±0.18 

1.64 ±0.09 

< (7 



Fig. 2.— SEDs of our faintest ALMA-selected galaxies — those detected in B6 {left) and B7 {right) - as a function of their possible 
redshifts. Their observed flux densities are indicated by the grey dots (un certainties smaller tha n the size of the dots). For this plot we 
have used the median SED for LABOCA sources from the LESS survey (ISwinbank et all 1201411 . Dashed lines sho w the limiting BAGS 
160-/rm and SPIRE 500-/im flux densities for one of the deepest survey carried out with Herschel (lElbaz et al.ll20nri . If they lie at z > 2, 
the faintest galaxies detected at the present stage of our survey would not have been detected in the deepest observations carried out with 
Herschel. They likely represent, therefore, a newly discovered population of galaxies, the bridge between Herschel- and SCUBA-detected 
galaxies and the classical UV-selected population. 


jets, though DSFGs at very high redshift, or very cold 
galaxies, may also have < Sbg due to the dust emis¬ 
sion peak shifting to redder wavelengths; we have not ap¬ 
plied such a cut. We have, however, discarded from our 
sample any submm detections with S'b 3 ,b 4 > 5'b6,B7- A 
clear example of the possible confusion between jets and 
DSFGs is ALMACAL J010838.75-b013455.9, one of two 
submm detections around the calibrator, J0108-I-0135. 
These two submm detections were first identified us¬ 
ing ALMA B6 data and they both resembled DSFGs. 
However, ALMA multi-band data revealed that AL¬ 
MACAL J010838.75-1-013455.9 is not detected in B7, 
contrary to the expectations for the SED of a high- 
redshift DSFG, despite the fact that B7 noise level would 
have allowed a detection. It also exhibits S'a.omm > 
<S'i. 2 mm, indicating that it is likely related to jet emis¬ 
sion from the calibrator. It is therefore excluded from 
our final sample of DSFGs in Table [21 On the other 
hand, ALMACAL J010838.56-b013504.2 is confirmed to 
be a distant DSFG due to its B7/B6 flux density ratio 


and its non-detection in B3. Two other clear examples of 
contamination lie around J1733—1304 and J0522—3627; 
despite lying far from the calibrators, high <S'b 3 , 4 /S'b 6,7 
ratios are clearly indicative of non-thermal emission. 

Turning now to a second caveat, over-densities of 
SMGs have been found around bright radio galaxies and 
quasars on a variety of y ales (e.g. Ilvison et al.l 1200(1 
l2008t IStevens et ^12003112004D . It is reasonable, then, 
to worry that our ALMA calibrators might also be re¬ 
lated to DSFG over-densities. In defense of our work, the 
properties of our calibrators - in particul ar their mass - 
are n ot as extreme as radio galaxies Isee lSevmour et HI 
120071) . which are the most massive galaxies at any red- 
shift, and which in any event are associated with over¬ 
densities which have proved difficult to prove conclu¬ 
sively; our calibrators are blazars, which are bright be¬ 
cause we are looking down the throat of their jets, not 
because they are hosted by spectacularly massive galax¬ 
ies. Nevertheless, we concede that selection of DSFGs 
around bright submm calibrators may be biased at a low 
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level. Indeed, a stronger bias is present in most if not 
all previous reports of faint ALMA counts in the liter¬ 
ature, since those previous works often derive number 
counts using ALMA science data that were deliberately 
centered on IR-bright or otherwise extreme galaxies that 
might also be associated with over-dense regions. 

An entirely unbiased determination of submm number 
counts (see lj4|) will only be possible when we can image 
large areas of the sky. While deep maps are possible 
with ALMA, wide observations are time-consuming due 
to the small FoV of the ALMA antennas in the bands 
that probe dust emission from star-forming galaxies. 

3.5. Spurious sources 

Creating large samples using a low SNR threshold, 
SNR ^3 — 4, leads inevitably to the inclusion of spu¬ 
rious sources. Working with samples contaminated by 
spurious sources can be a valid approach if our interest 
lies only in the study of number counts, since we can 
determine and apply accurate correction factors. 

Here, instead, we have opted to be conservative - 
to reduce as much as possible the contamination from 
spurious sources - since one of our aims is to study 
the nature of the detected DSFGs. For this we need 
to be sure that the sources we detect are real, lest 
we were to identify a spurious source as a dropout 
in another band, perhaps initiating a long and sorry 
saga. In our source selection, therefore, we have se¬ 
lected sources at SNR > 5, a threshold higher than has 
been adopt e d in a l l recent reports on faint ALMA counts 
llOno et all I20l4 iFuiiinoto et al.l l20I5t iCarniani et al.l 
I20I5I : iHatsukade et al.l 1201311 . While this will signifi¬ 
cantly decrease the contamination from false detections 
in our final sample of DSFGs, we have inevitably ex¬ 
cluded some of the faintest DSFGs from our catalogue. 
While we can be reasonably confident that the number 
of false detections in our resulting 5cr sample should be 
zero, we can state with absolute confidence that all six 
of the DSFGs with multi-band detections, as discussed 
above and shown in Table [U are real. 

3.6. Completeness 

We have employed the traditional method of injecting 
simulated sources with different flux densities at random 
positions of our ALMA maps in order to study the com¬ 
pleteness of our selection method as a function of SNR. 
All previous works in this area have added simulated 
sources to the clean maps. Since we are working with 
interferometric data we have instead opted to perform 
our simulations in the visibility plane. We have taken 
ALMA calibrator maps where no DSFG is detected and 
inserted artificial, point-like sources with flux densities 
ranging from O.I to lOmJy in steps of 0.025 mJy, thus 
covering a wide range of SNRs. Sources were always 
injected within 1.5 x fwhMpb, the same area used for 
source detection (E3D- We repeated the source injec¬ 
tion 100 times for each value of the input flux density. 
Once a source has been injected at a given position, we 
shifted the phase center of the map to those coordinates 
to ensure the source should show up at the center of 
the dirty map. Because this is done for all the injected 
sources, it is possible to define exactly the same cleaning 
box for all of them, ensuring that the cleaning, source 



Fig. 3.— Ratio between flux density at 870/.tm and 1.2 mm for 
the six DSFGs with ALMA detections in both B6 and B7. The 
shaded area represents the ratio expected for the composite SED 
of SMGs from the LE SS survey for the redshift range, z = 1-3 
(ISwinbank et al.l 120141) . High S 870 /im/'Si. 2 mm ratios can be ex¬ 
plained with a higher dust emissivity index. Low S 870 fjm/Sl. 2 mm 
ratios could be indicative of higher redshifts. 

detection, and flux-density determination is a homoge¬ 
nous and consistent process. Imaging and cleaning was 
accomplished in non-interactive mode, down to the r.m.s. 
of the dirty map. As will be shown in H3.81 this clean¬ 
ing method provides accurate values of the recovered flux 
density of the simulated sources and also represents the 
optimal method for measuring the flux density of our de¬ 
tected DSFGs. We consider than an injected source has 
been recovered if it is detected with SExtractor at > 5a 
within a synthesised beam of the center of each map. 
Glean maps, prior to PB correction, are used for source 
detection, consistent with the process used to select our 
DSFGs (see ^3.1|) . 

As a result of our simulations we have determined 
that our survey is nearly 100% complete at SNR > 7 
and ^80% complete at SNR > 6. These SNR values 
are higher than those reported in most previous works 
that have analysed faint DSFGs detected in deep ALMA 
maps. For example, the completeness in I Simpson et al.l 
(|2015cll is 93% at > 4g, r i sing t o about 100% at 5.5cr, 
meanwhile iFuiimoto et ahl (|2015l l claim their ALMA ob¬ 
servations are 90% complete at SNR > 4.5 (their Fig¬ 
ure 3). The higher completeness at a given SNR de¬ 
rived in previous works is a consequence of the lower 
SNR threshold used for source selection. Previous works 
have employed selection thresholds at SNR > 4, or lower, 
which increases the completeness at low flux densities (or 
SNRs) at the cost of increasing the number of spurious 
detections. 

3.7. Effective area 

The sensitivity of any single-pointing ALMA image (or 
any interferometric image) decreases with increasing dis¬ 
tance from the center of the map due to the PB response 
of the individual antennas. The effective area sensitive 
to a given flux density therefore varies with flux density. 
As an example, a galaxy detected in the center of the 
map at 6a would not be detected at > 5a were it located 
instead at the edge of the map (defined by 1.5 x fwhm) 
and would not then be included in our sample. 

Obtaining the relationship between the effective area 
of our survey and the flux density of the sources detected 
is important since we have to correct the number counts 




























Area [arcmin 


ALMACAL I: a deep and wide submm survey 


9 



Fig. 4.— Effective area covered by our ALMA submm survey as 
a function of the flux density of the detected sources for the two 
bands considered in this work. The sensitivity of a given ALMA 
calibrator map is a function of the distance to the center of the 
map. The effective area covered by the survey is then a function 
of the flux density of the detected sources. For a given source 
detected at a given band with a given flux density, we calculate 
the area where that galaxy coul d ha ve been detected at > 5 (j (our 
limit for source detection - see ED in all pointings in that band. 
In this calculation we have modelled the PB response (given by 
CASA during imaging) with a Gaussian. 

for this effect (see Q. For a given flux density, Sin, 
and a map with an r.m.s. of Cin we calculated the SNR 
at the center of the map, SNRcen = Sin/cin- We then 
obtained the radius, rum, at which that SNR decreases 
to five, which gives the area in which a galaxy with a flux 
density, ^in, can be detected. If 2 x rum exceeds 1.5 x 
FWHM, we define rum from 2 x Him = 1.5x fwhm. We 
performed this calculation for all observations in B6 and 
B7 independently, since we will report number counts 
for both bands. Figure |4] shows the effective area of our 
survey as a function of the flux density of the detected 
sources in B6 and B7. The total number of pointings in 
B7 is smaller than in B6; taken together with the smaller 
area of each B7 observation, this means that the total 
area covered in B6 (around 16 arcmin^) is considerably 
higher than that covered in B7 (around 6 arcmin^). 

3.8. Flux boosting 

It has long been known that the flux densities of galax¬ 
ies detected at relatively low SNR can be boosted due to 
the presence of noise fluctuations. We have analysed the 
effect of flux boosting for our ALMA detections using 
the same set of simulated point-like sources used in 113.61 
Once the phase center of the map has been shifted to 
the position of each injected source, and the visibilities 
have been imaged, as explained in 113.61 the flux densities 
of the detected galaxies are measured in the clean maps 
using IMFIT prior to correction for PB attenuation. 

We find that the flux density is over-estimated by up to 
20% for sources detected at SN R ^ 5, which is in goo d 
agreement with the findings of iSimnson et al.l (l2015cf) . 
By contrast, o ur results are in marke d disagreement with 
the findings of iFuiimoto et al.l (|2015ll . where 5'out/5'in ra¬ 
tios are reported to be almost unity with a variance of 
only 5%. They argue that their low contamination from 
flux boosting is due to the high spatial resolution of their 
ALMA im ages. However, we rea ch a similar spatial res¬ 
olution as iFuiimoto et al.l (|2015D and find that the ef¬ 
fect of flux boosting i s sign ificant at low SNR. Addition¬ 
ally, IFuiimoto et al.l (|2015ll considered detections down 


TABLE 3 

Cumulative number counts derived in our work. 


S'sTOum [mJy] 

-/V(> 5'870Mm)[xl0^deg 

0.4 


1.0 

3-8^:? 

'S'l.2mm [n^Jy] 

N{> Sl.2mm)[xl0^ deg-^] 

0.2 

5-61I 

0.8 

0^ 


to SNR ~ 3.7, where the effect of the flux boosting should 
be >30% if we extrapolate from our > 5a results to lower 
SNR. Many of our ALMA-detected DSFGs were detected 
in the SNR regime where the effect of flux boosting is not 
significant (see Table [2|). The flux densities of our DS¬ 
FGs detected at SNR < 10 have been corrected for the 
boosting effect. 

The analysis of the ratio between input and output 
flux densities of the injected sources also shows that 
‘S'out/'S'in ~ 1 for sources detected at > lOcr, confirm¬ 
ing that our non-interactive imaging method provides 
robust flux densities. On the other hand, we obtain that 
the variation in S'out/'S'in from galaxy to galaxy at low 
SNR can be significant (up to a factor of two of differ¬ 
ence), meaning that caution should be exercised when 
interpreting individual flux density ratios. 

4. NUMBER COUNTS 

In this section we present the cumulative number 
counts obtained for our survey of DSFGs and compare 
them to previous results published in the literature. 
Since we have a reasonable number of submm detections 
in B6 and B7, we present for the first time ALMA cu¬ 
mulative number counts for each band independently. 

The contribution to the cumulative number counts of 
a source, i, with a flux density. Si, is: 


N,{S^) 


1 /sp(Si) 


( 1 ) 


where /sp(<S'i) is the fraction of spurious sources at Si, 
C{Si) is the completeness of the survey at Si, and A{Si) 
is the effective area covered by our survey at Si. The 
completeness and effective areas at different flux densi¬ 
ties are taken from 113.61 and respectively. Thanks 
to our multi-band and multi-epoch observations and our 
conservative selection criterion, the fraction of spurious 
sources is expected to be zero: fsp = 0 for any Si. In 
order to calculate the cumulative number counts at dif¬ 
ferent flux densities, we have to sum over all the galaxies 
with flux densities higher than the adopted values: 


N{>S)= 

ySi>s 


1 - /sp(g0 

G(5,) ■ A{S,) 


( 2 ) 


Figure [S] shows the cumulative number counts de¬ 
rived in this work at 1.2 mm and 870/rm, where the 
flux boosting, completeness, and effective area fac¬ 
tors have been included (see also Table o. For each 
band, we have reported number counts at flux densi¬ 
ties that evenly divide the sample into two bins with 


similar numbers of galaxies. Previous re sults based on 
ALMA data, from iKarim et alJ (l20I3fl. iSimpson et al.l 
(|20I5cll . iHatsukade et all ( 201311 . lOno et aLn20I4ll and 
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Fig. 5.— Cumulative number c o unts of DSFGs der i ved in this work a t 870/xm (left) and 1 . 2 mm {right). We also plot the results 
reported by IHatsukade et all 1120131) ; IKarim et al.l 1120131') ; lUno et II2014I) ; IFujiruoto et al.l II2015I) ; ISimpson et al.l ll2015d) . Number counts 
determined in previous work at different wave lengths have been converted into 870- /rm and 1.2-mm count s assuming the typical SED of 
bright SMGs at ^ = 2.3 HSwinbank et al1l2014l) . The grey curve is the fit obtained in ISimpson et ahl ll2015d) extrapolated towards the fiux 
densities covered in our work . Our counts are lowe r than those presented in previous works; spe cifically, our counts are lower by a factor 
of around 2x with respect to lFuiimoto et ^ II2015I) and by a factor of about 7x with respect to IHatsukade et~a n <2013) . 


iFuiimoto et all (|2015D are also shown. We have decided 
not to compare with cumulative number counts derived 
with ground-based single-dish telescopes or Herschel 
since bright SMGs are sometimes resolved into multiple 
components, which may strongly affect the deriyations 
of nu mber counts (jKarim et al.1 1201, 'll : iSimpson et al.l 

The number counts deriyed in previous work have 
been determined in different ALMA bands. Cumula¬ 
tive number counts derived from ALMA follow-up ob¬ 
servations of SCUBA- or LABOCA-selected SMGs have 
been determined in B7. Meanwhile, previous work re¬ 
porting faint ALMA counts using archival ALMA sci¬ 
ence data have largely been focused on B6 data due 
to the larger instantaneous FoV. In order to avoid un¬ 
certainties due to converting flux densities between dif¬ 
ferent bands, we have compared only results obtained 
in the same ALMA bands, with one exception. In or¬ 
der to include bright counts in B6. we also include 
the re sults from iSimpson et all (|2015cll and IKarim et al.1 
(|2013f) . The counts have been converted from 870/rm 
to 1.2 mm us ing the composite SM C SED from the 
LESS survey (jSwinbank et al.l [20141) . redshifted to z = 
2.3, the median red shift of LABOCA-detected SMGs 
(jSimpson et al.ll2014f). This conversio n factor is the same 
as that used iniSimpsoii et al.l (|2015cll and similar to that 
used in lEuiimoto et al.1 (l2015f). We have simi l arly c on- 
verted the 1.3-mm number in IHatsukade et ^ (l2013f ) to 
1.2 mm. 

Any comparison of number counts between different 
works will, of cour se, be affected b y thi s transforma¬ 
tion. The SED in iSwinbank et al.1 ( 20141) was derived 
for bright SMGs, whereas many of our DSFGs are faint 
and might thus have different dust temperature distri¬ 
butions. Also, the median redshift of the faint DSFGs 
in our sample could differ from the bright sample of 
iSimpson et al.l (j2014D . and the redshift distribution also 
likely de pends on the wavel ength where the DSFGs were 
selected (iVieira et al.ll2013l) . 


It can be seen from the left panel of Figure [S] that 
the cumulative 1.2-mm number counts we obtain at 
the present stage of the survey are lower than previ- 
ous claims. Ev e n for the case with the best agreement, 
IFuiimoto et al.1 (j201,5D . we obtain lower number counts 
by a factor of 2x. A more s evere disagreement is found 
with IHatsukade et al.l (I2013D , who report number counts 
significantly higher than ours. 

One of the main differences between our work and pre¬ 
vious work is the source-detection procedure. We se¬ 
lected only galaxies detected at > 5a, while all previ¬ 
ous work reporting faint number counts have included 
sources detected at lower SNR. This might lead to the 
inclusion of spurious detections, as each of these pre¬ 
vious papers pointed out, whereas our process leads to 
negligable contamination. A fair comparison with pub¬ 
lished number counts would require consideration of only 
the > 5cr detections in those works. However, most do 
not give the flux density, SNR, area covered, and com¬ 
pleteness for individual detections. Furthermore, com¬ 
parison with previous work requires a re-calculation of 
the area covered and completeness of previous studies, 
since th ese depend upon the source detection criteria em¬ 
ployed. lOno et ahl (j2ni4[ ) did report the list of detections, 
including their SNR. Among their 11 submm detections 
in ten ALMA maps, only two are detected at > 5cr, only 
one of them as bright as the DSFGs in our sample. 

The right panel of Figure [5] represents the cumulative 
number counts derived in this work at 870 /rm. Since all 
previous work studying faint number counts with ALMA 
has been focused on lower frequencies to take advantage 
of the larger instantaneous FoV, our results represent 
the first determination of faint number counts at 870 /rm 
with ALMA. I n Figure [5] we combine our counts with 
those d erived inlS imoson et al.l ()201,5cD and IKarim et al.1 
(1201,31) . ISimpson et al.1 (|2ni,^cD fi tted a double power law 
to their own data, and those of IKarim et al.1 (|2ni,3D . It 
can be seen that the extrapolation of their best-fit func¬ 
tion toward S'sTO/xm ~ 0.3 mJy reproduces the cumulative 
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number counts obtained in this work. 

Our wide and deep submm survey, exploiting ALMA 
calibrator data, is proving to be very powerful for deriv¬ 
ing faint number counts, both in ALMA B6 and B7. This 
is thanks to the low noise levels that can be reached when 
combining data for a specific calibrator and the area that 
can be covered by combining deep data for different cal¬ 
ibrators. As more data become available for a larger 
number of calibrators, we expect to derive more robust 
number counts in B6 and B7, based on larger samples, 
obtained with deeper observations over larger areas. The 
large number of calibrators used during ALMA observa¬ 
tions will enable us to determine robust number counts 
even in B8, despite the relatively small FoV. The combi¬ 
nation of number counts in all these bands will provide 
strong constraints on models of galaxy formation and 
evolution. 

5. CONCLUSIONS 

In this paper we have presented a novel technique 
which exploits publicly available ALMA calibration data 
to carry out a deep and wide submm survey, sensitive to 
dusty, star-forming galaxies. Our main conclusions are: 

1. We have demonstrated the power of using ALMA 

calibrators to detect and study high-redshift DS- 
FGs by analysing relatively deep data, down to 
r.m.s. ^ 25 /iJy beam“^, of 69 calibrators. We have 
focused on ALMA bands 6 1.2 mm) and 7 

(~ 870 /rm), since they probe the dust emission of 
high-redshift star-forming galaxies. We have cov¬ 
ered about 16 and 6 arcmin^ in ALMA B6 and B7, 
respectively. 

2. Using a conservative detection threshold, > 5cr, 
we have found a sample of eight and 11 DSFGs 
in ALMA B6 and B7, respectively. Among these, 
six are detected in both B6 and B7 and the flux 
density ratios between those bands are consistent 
with high-redshift DSFGs. Six of our DSFGs have 
ALMA B3 (3 mm) observations, where they are un¬ 
detected, again compatible with high-redshift DS¬ 
FGs. 

3. The average 870-/rm flux density of our DSFGs 
is lower than those of the classical population of 
single-dish-selected SMGs, and comparable with 
the (stacked) flux density of LBGs at z ^ 3. 
The faintest galaxies detected in our survey would 
have been missed even by the deepest Herschel 
extragalactic surveys, for any reasonable redshift, 
demonstrating the relevance of our survey for 
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studying the faint population and for exploring the 
link between the extreme population of SMGs or 
Herschel-selected galaxies and the more abundant 
population of UV-selected galaxies. 

4. Using our sample of DSFGs we have determined 
cumulative submm number counts at 870 ^m and 
1.2 mm, independently, from a sparse sampling of 
the astronomical sky, thus remaining relatively free 
of cosmic variance. This is the first determination 
of faint number counts with ALMA at 870 /im. We 
find that the counts are lower than previously re¬ 
ported by a factor of at least 2x at the lowest flux 
densities probed by our survey. 

We plan to increase the depth of the survey and the 
size of the area covered, as observations of more calibra¬ 
tors become available. Besides the obvious advantage of 
exploiting the significant fraction of ALMA time that is 
required for calibration, our approach has several other 
key advantages. Multi-epoch, multi-band observations 
mean that our detections of DSFGs are entirely secure, 
with confirmation via appropriate submm spectral in¬ 
dices, at matched resolution, and that their redshifts will 
ultimately be determined via the blind detection of spec¬ 
tral lines. Uniquely, our approach enables morphological 
studies of faint DSFGs - a more representative popula¬ 
tion of star-forming galaxies than conventional SMGs - 
in fields where self-calibration is feasible on timescales of 
a few seconds, meaning that baselines yielding milliarc- 
second spatial resolution can be exploited, for free. 
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